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Abstract 

We find the condition that the primordial baryon asymmetry is not washed out in the SU(2)l 
triplet Higgs model by solving the Boltzmann equation. We further require that the model is com- 
patible with the recent results of the neutrino oscillation experiments and the Wilkinson Microwave 
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obtain the allowed region of the parameters in the model. 



Anisotropy Probe, and the constraints on the p parameter imposed by the CERN LEP. We finally 
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I. INTRODUCTION 



The executed experiments of the neutrino oscillation have made it certain that the neutri- 
nos have the very small, but finite masses 

mm 

. The recent KamLAND experiment 
gave us the more precise information about the neutrino mass matrix 

Q,y,3- Although 

the neutrino oscillation experiments are not able to determine the magnitudes themselves of 
the neutrino masses, the recent results of Wilkinson Microwave Anisotropy Probe (WMAP) 
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Since 



imposed the upper bound of the sum of the neutrino masses 
the neutrinos have no masses in the standard model, it is an urgent issue that we extend 
the standard model and build the model where the neutrinos naturally have the very small 
masses in comparison with the other leptons and quarks. It is expected that the smallness 
of the neutrino masses is explained by the Majorana type mass, which only neutrinos among 
quarks and leptons are able to have. The Majorana mass terms have the noticeable feature 
that they violate the lepton number, which is conserved at the tree level in the standard 
model. A few types of models of the Majorana neutrinos are proposed. The seesaw model 
3, 16, Q, the SU(2) L triplet Higgs model and the Zee model Q are the three 

typical models where the standard model is minimally extended to induce the Majorana 
masses of the neutrinos within the framework of the SU(2)l ® U(l)y gauge theory. 

In this paper we focus on the SU(2)l triplet Higgs model. Especially, we introduce 
the model which is built in the framework of the large extra dimension and can naturally 
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explain the smallness of the neutrino masses by the "shining" mechanism 
in one brane and the lepton number is violated in the distant other brane. A messenger 
field informs us of the lepton number violation in the other brane. Since the messenger field 
spreads out in the extra dimensions, the density of the field is diluted, which is called the 
" shining" mechanism. Since we are informed of the only tiny lepton number violation, the 
very small neutrino masses are naturally induced in our brane. Furthermore, based on the 
't Hooft's naturalness condition, the tiny lepton number violation, that is, the very small 
neutrino mass, is stabilized under the quantum corrections. In this paper we adopt the 
4-dimensional effective theory in our brane of the SU(2)l triplet Higgs model in the large 
extra dimension. It is also noticed that the SU(2)l triplet Higgs model must not induce 
the dangerous pseudo Majoron which is inconsistent with the accelerator experiments as the 
CERN LEP, and the model is constrained from the observation of the p parameter in the 
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LEP [23, 24]. The lepton flavor violation in this model is analyzed in |25l l26j]. 
The baryon number in the present universe is well known to be [10| 

B = ^ = (8.4 ~ 9.1) x 10~ n , (1) 

where ns is the baryon number density and s is the entropy density. Some scenarios where 
the baryon number is dynamically gen erated from vanishing initial baryon number, Bi n i = 0, 
are proposed. The GUT scenario |27j,^a, the electroweak scenario |29j, lJQJ], the Fukugita- 

n 

Yanagida scenario [31], and the Affleck-Dine scenario [22J are proposed as the typical sce- 
narios of the baryogenesis. The GUT and electroweak scenarios have already been almost 
ruled out. It is often mentioned that these scenarios of the baryogenesis have a double-edged 
behaver, that is, the scenarios of the baryogenesis simultaneously have a potential to erase 
the existing baryon number. Especially, while the mass models of the Majorana neutrinos 
have a potential to generate the baryon number through the leptogenesis as in the Fukugita- 
Yanagida scenario, these models simultaneously have a potential to erase arbitrary initial 
lepton and therefore baryon numbers. There exist the sphaleron-anomaly processes which 
create and annihilate the baryon in an equal amount as that of the lepton in the standard 



model 
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371 ] . When the the sphaleron-anomaly processes and the lepton num- 



ber violating processes in the mass models of the Majorana neutrinos are simultaneously 
in equilibrium, the arbitrary initial baryon and lepton numbers are actually washed out. 
Since the sphaleron-anomaly processes are in equilibrium at the large region of tempera- 
ture between 100 GeV and 10 12 GeV, whether the lepton number violating processes are in 
equilibrium or not in the region handles the fate of the cosmological baryon number. 

In this paper, we assume that the primordial baryon number is generated by some un- 
known mechanism at the very early universe and require that the primordial baryon number 
is not washed out in the mass models of the Majorana neutrinos where the baryon number 
is never generated. The situation where the baryon number is never generated via lepto- 
genesis in the Majorana mass models is typically realized when we assume that the lepton 
sector has no CP violation because of the Sakharov's three conditions [38j. Although we 
assume no CP violation in the following analyses, the results of this paper are valid as long 
as the mass models of the Majorana neutrinos are not able to generate the present baryon 
number. Actually, it is shown that even if the SU(2)l triplet Higgs model, which we use in 
this paper, has CP violation in the lepton sector, the baryon number can not be generated 



3 



40 



4l|. 



3^1 . The leptogenesis in the other extended triplet Higgs models is analyzed in 
In order to generate the baryon number via leptogenesis, two triplet Higgs fields must be 
introduced. The similar analyses are carried out in the GUT scenario |42l . and in the 



Fukugita-Yanagida scenario (44 , 
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In order to escape from washing out of the baryon number in the mass models of the 
Majorana neutrinos, it is necessary that the lepton number violating processes are out of 
equilibrium. An out-of-equilibrium process means that the reaction rate of a process is 
smaller than the Hubble parameter at that time. In order to violate the lepton number 
explicitly in the SU(2)l triplet Higgs model, two kinds of interactions must be newly in- 
troduced. At first sight, it is not clear whether both of the two kinds of interactions must 
be out of equilibrium or either of the two kinds of interactions must be out of equilibrium 
in order to protect the primordial baryon number. The answer we obtain in this paper is 
that if either of the two kinds of interactions is out of equilibrium, the primordial baryon 
number is not washed out. If either of the two kinds of interactions is out of equilibrium, 
the lepton number is considered not to be explicitly broken and an approximately preserved 
global U(l) charge which contains the lepton number exists. Thus the primordial baryon 
number is protected in proportion to the initial value of the approximate U(l) charge. It is 
the first purpose in this paper that we can concretely estimate the condition to protect the 
primordial baryon number by solving the Boltzmann equation in the SU(2)l triplet Higgs 
model. We next analyze whether the SU(2)l triplet Higgs model can simultaneously satisfy 
the obtained condition to protect baryon number and the results of the neutrino oscillation 
experiments or not. We further require that the model satisfies the results of WMAP and 
the constraints on the p parameter. It is the second purpose in this paper that we require 
that the model satisfies all the above conditions and obtain the allowed region of the pa- 
rameters in the model. The similar analyses are carried out in the seesaw model and in 
the Zee model J3, 1 . 

The outline of this paper is as follows: in Sec. HH we review the SU{2)l triplet Higgs 
model where the smallness of the neutrino masses is naturally explained. We make sure 
that the dangerous pseudo Majoron which is inconsistent with the accelerator experiments 
as LEP does not appear in the model. In Sec. Hill we obtain the condition to protect 
the primordial baryon number by solving the Boltzmann equation in the SU(2)l triplet 
Higgs model. In Sec. IIV| we require that the SU(2)l triplet Higgs model satisfies the 
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condition to protect the baryon number, the results of the neutrino oscillation experiments 
and WMAP, and the constraints on the p parameter. We finally obtain the allowed region 
of the parameters in the model. Sec. 0is devoted to a summary. 



II. SU(2) L TRIPLET HIGGS MODEL 

In this section we first introduce the SU(2)l triplet Higgs model in the normal four 
dimensional space-time as an effective theory of the model in the large extra dimension 
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It should be noted that two kinds of interactions must be newly introduced in 
order to violate the lepton number explicitly. We second show that the very small neutrino 
masses are naturally induced and simultaneously confirm that the dangerous pseudo Majoron 
which is inconsistent with the accelerator experiments as LEP does not appear under the 
suitably assumed Higgs potential. We finally check how the observations of the p parameter 
at LEP constrain the model. 

We newly introduce the SU(2)l triplet Higgs fields A assigned with the hypercharge 
Y = 2 in addition to the fields of the standard model as 

aJ (+ '^ ?++ V (2, 

\ e -e/V2j 

The Yukawa interactions of the triplet Higgs with the leptons are written as 

c yuka W a = -I/^TVp^^A] + h.C. (3) 



where 7}*^ is a SU{2) L triplet which is composed of the two lepton doublets, 



+ h.c. 



and f af3 is the symmetric Yukawa coupling constants (f af3 = f 13 "). Since we take the base 
where the mass matrix of the charged leptons is diagonalized, the indices, a and (3 take e, p, 
or r. When the neutral triplet Higgs field £° gets the vacuum expectation value (VEV), 
the interactions in Q, whose mass dimension is four, induce the Majorana mass terms of 
the neutrinos, which explicitly violate the lepton number. However when we assign the 
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triplet Higgs fields with the lepton number La = +2, none of the interactions in the model 
violate the lepton number at the symmetric phase where the neutral triplet Higgs field 
has no VEV, that is, the lepton number is not explicitly violated. In the broken phase, 
neutrinos acquire Majorana masses via the spontaneous breakdown of the lepton number, 
which in turn means that a Majoron, which is inconsistent with the accelerator experiments 
as LEP, appears. Thus, in order to avoid the dangerous Majoron, we need to violate the 
lepton number explicitly. For such purpose, we further introduce the cubic interaction of 
the triplet Higgs with ordinary Higgs doublet $ = (0 + , 0°)*, 

C cubic = ~ATr[T*, 9 Ai] + h.c. 



-'-A 
2 



^) 2 r"-v / 2rt o r-(0°) 2 rj+h-c., (5) 

where 7$,$ is a SU(2)l triplet which is composed of the product of the Higgs doublet : 



, I "* V ^ I , (6) 

* -0°0° + 0" 

which has the hypercharge Y = +2 and the lepton number L = 0, and A is the coupling 
constant. 

The coupling constant A in (|SJ) can be naturally made very small in the large extra dimen- 
sion scenario by the shining mechanism j^ll . I22I ] . Since the neutrino masses are in proportion 
to the VEV of the neutral triplet Higgs field and the VEV turns out to be in proportion to 
the coupling constant A, the smallness of the neutrino masses is naturally explained by the 
smallness of the coupling constant A in the model. Here, under a simplified Higgs poten- 
tial, we show that the VEV of the neutral triplet Higgs field is in fact in proportion to the 
coupling constant A and simultaneously confirm that the dangerous pseudo Majoron does 
not appear. We work in a simplified Higgs potential, 

V($, A) = -// 2 $t$ + A($t$) 2 + M 2 Tr[A t A] + -(A Tr[T^A<) + h.c), (7) 



where /1 2 and M 2 are positive . We extract the sector with electromagnetic charge Q = 
in 0, 

u 2 A M 2 

V(<Pu 2 , 6, 6)q=o = - y(0? + <Pl) + ^(01 + <Pl) 2 + — + 

A +20i0 2 6l, (8) 



2V2 



where the real fields (0i, 02, £1,^2) are defined as 



,0 _ 0i + #2 , t0 _ + 



Solving the conditions for the extremum of the potential (|5J). 

^Vq=o = dVQ=o = <9^q=o = <9Vq=o = 

(901 <90 2 <9£i <9£ 2 

we obtain the vacuum expectation values of 0i, 02, £1, and £2 



and C = S1 ' s " . (9) 



0. (10) 



u Av 2 

< 0i >= , < 02 >= 0, < 6 >= — -= , and < £ 2 >= 0, (11) 

V v/A-A 2 /4M 2 ' ^ U 2v/2M 2 ' W 1 ; 

and define as t> =< 0i > and i>a =< £1 >. The shifts, a, N, a A , Aa from the VEVs (fTTj) are 
defined as 



n, V2 



{v + a, N) and &) = (v A + a A , N A ). (12) 



Accordingly, the potential (jBJ) is rewritten as 

V(a,N,a A ,N A ) Q ^ = -\v' - ^ + i(2A„V + " 

+ ^ + ^ ~ <13) 
where we have kept only the terms up to quadratic in the fields. Here we define the potential 
V{a,a A ) and V(N, N A ) as 

V(a,a A ) = ^(2Xv 2 )a 2 + ~M 2 a 2 A - 4=<*t a , 

V(N, N A ) = ^N 2 + \M 2 Nl - ^=NN A . 

The potential V(a, a A ) is diagonalized by an orthogonal transformation 

a I = R{0a) I ° I , (14) 

where the rotation matrix R(8 a ) is defined as 

. cos 9 a — sin o . 
R{9„) ={ I ■ (15) 

sin cos v a 
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Here the rotational angle 8 a and the two mass eigenvalues, m 2 , and m 2 , are given as 

^ - «,.-Cv (16) 

„ 9 M 2 + 2Xv 2 ± v/(M 2 - 2Xv 2 ) 2 + 16M 2 v 2 A /v 2 . , 

w CT /,m^ = , (17) 

where m 2 , takes the positive sign. In a similar way, in order to diagonalize the potential 
V(N, we rotate the base (N, N A ) as 

N' 
N' A 

where the rotational angle 6 N and the two mass eigenvalues m? N , and m 2 N , are given as 




(18) 



^"■v = ^T^' ( 19 ) 



v 2 — Av A 



m^, = and m^/ = M 2 ( 1 + 4-^ ) . (20) 



We next extract the terms with non-zero electromagnetic charge, Q ^ in (JJJ) : 



Vq # o = Vl0 + | 2 + A(|0 + | 4 + 2|0+| 2 |/| 2 ) + M 2 (\t? + K ++ | 2 ) 

"^((0 + ) 2 r _ - v^vr ) + h-cl . (2i) 



i 

+ 2 

Since the potential (}2Tf induces the mixing term between <f) + and £ + at the symmetry broken 
phase, we diagonalize the mixing term by using a rotation, 

= R{B+) K l) ■ (22) 
Here the rotational angle 9 + and the two mass eigenvalues Tn^ >/+ and m 2 l+ are defined as 

tan2£ + = -^jr, (23) 




m 2 ,^ = and m 2 , + = m 2 , = M 2 (l + 4%^) . (24) 

From (f2T)|) and (|2*4^1 we thus confirm that there are no dangerous pseudo Majoron which is 
inconsistent with the experiments. In fact, N' and (ft + are would-be N-G bosons, and the 
masses of remaining physical Higgs scalar N' A and £' + can be kept large in the case that M 2 
is large enough, even if A is small. When the gauge symmetry is spontaneously broken and 



the neutral Higgs field, £° gets the VEV VaJ V%, the Majorana mass terms of the left-handed 
neutrinos are induced in (j3J) as 

C Majorana = --^(^I/j? + h.C., (25) 

where the neutrino mass matrix m al3 is defined as 

m* = r < e >= r#^. w 

We naturally expect that the magnitudes of the coupling constants f a/3 are between the 
Yukawa coupling constants y e and y T of e and r, that is 

y e ~ 2.9 x 10- 6 < |r^| < y T ~ 1.0 x lO" 2 . (27) 

The coupling constant A is expected to be as small as the order between 1 eV and 1 MeV 
in the models of the large extra dimension at classical level 3, The masses of the 
triplet Higgs fields, M is expected to be of the order of the electroweak scale. When we fix 
A and M at A — 1 keV and M — 1 TeV, for example, the neutrino masses are of the order, 
10~ 5 ~ 10 _1 eV, which can be enough small to satisfy the results of the neutrino oscillation 
experiments and WMAP. Furthermore, since the lepton number is conserved in the limit 
A — > 0, the naturalness proposed by 't Hooft guarantees the stability of the small A under 
the radiative corrections. Thus, we make sure that the smallness of the neutrino is naturally 
explained in this model. 

We finally check how the observations of the p parameter restrict the 577(2) triplet Higgs 
model. The p parameter is defined and calculated to be 

M| _ , 2 + 2,| vj 
f M 2 z cos 2 6 w v 2 + Av\~ v 2 ' { ' 

where 9w is the Weinberg angle and Mw and Mz are the masses o: 



the W and Z gauge 
23l 2^1 restricts the ratio 



bosons. The constraint on the Ap = p — 1 by LEP experiments 
of VEVs, va/v, as 

- = -4— < 0.03. (29) 
v 2V2M 2 V ; 

It is noted that the constraint (|29|) gives the upper bounds on the coupling constant A for 
the fixed M, but does not restrict the coupling constant f a/3 at all. 
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III. WASHING OUT OF BARYON NUMBER 



We want to figure out the condition that the initial baryon number is not washed out in 
the model we consider, assuming that the initial baryon or lepton number was generated by 
some unknown mechanism in the very early universe. In this section, we obtain the condition 
in the SU(2)l triplet Higgs model by solving the Boltzmann equation. In subsection IIII Al 
we construct the Boltzmann equation in the SU (2) l triplet Higgs model for the case of no 
CP violation in the lepton sector. In subsection IIII Bl we solve the Boltzmann equation 
and derive the time evolution of the lepton and baryon numbers. In subsection IIII CI we 
obtain the condition for the baryon number not to be washed out. The way in which the 
baryon number is washed out in the minimal SU(5) GUT with no CP violation is analyzed 



in 



42, 



43|. 



A. Boltzmann equation in the SU(2)l triplet Higgs model 

Since our purpose is to know the time evolution of the difference between the particle 
number and the anti-particle number, we define the baryon number B, the lepton number 
L, the charged gauge boson number W, the doublet Higgs number $, and the triplet Higgs 
number A as 



B EE — 



W 
$ 
A 



Hb n b -n- b 
s 



E£ \ n u iL n u iL n u iR n u iR 
31 S S 

t=u,c,t 



i=d,s,b 



1 \ n dtL - n- A%h | n dm - n diR 
3 L s s 



ul _ ni - nj 
s s 



E 



a=e,fi,r 



n e a L n e a L _|_ n e a R n e aR 



n w _ n w - - n w + 

1 

s s 



+ 



(30) 



(31) 
(32) 
(33) 
(34) 



s s s s 

where nx is the number density of a particle X and s is the entropy density in the universe. 
We rewrite the numbers defined above by using the relation described by chemical potential 
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U) as 



i=u,c,t i: 

L = k x (^„ a + + 



i=d.s.b 



W = k x 4/i w , 



(35) 

(36) 
(37) 



$ = k x 2(-^- + /ty>), (38) 
A = k x 2(-/i ? — - + /i^o), (39) 

where /ix is the chemical potential of a particle X and k is defined as k = T 2 /3s. We next 
define the total electromagnetic charge Q and the total third component I 3 of the SU(2) l 
isospin as 

rii — nj \ \ -» 1 \ -» 

i = k x |^ + 2 Vi/^ifl'j 

i=all i=boson i=fermion 

k x - 4//£— - 2/^- - 4/^- 



Q = ^Qi 



(40) 



i=all 



i=u,c,t 

■ rij — m 



i=d,s,b 



k x 



"/V 



/i^O — 2/l£-- — 2/lgO — 4yU, t 



(41) 



i=u,c,t i=d,s,b a=e,[i,T 

where and Jg are the electromagnetic charge and the third component of the isospin 
assigned to the particle i, and gi is the internal degree of freedom of the particle. 

The time evolution of the particle numbers in a model, of course, is caused by the inter- 
actions in the model. We enumerate each interaction in the model below. The interactions 
(i) ~ (vi) listed below exist in the standard model. 

(i) The gauge interaction of the quarks : C = g / y^un^V^djLW^ + h.c. 

Here Vkm is the Kobayashi-Maskawa matrix. This Lagrangian density causes the 
processes, 



d 



UiL + W-. 



(42) 



(ii) The gauge interaction of the leptons : C = g/ ' sf^VoiLl^^aiWt + n - c - 



e aL < — ► v aL + . 



(43) 
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(iv) The Yukawa interactions between the up-type- quarks (down-type-quarks) and the dou- 
blet Higgs fields : C = y l d d iL d iR (j) + y l u u iL u iR (f) * + h.c. 



=> d iR + 0° < > diL, u iR < — > u iL + 0°. 

(v) The Yukawa interaction between the leptons and the doublet Higgs fields : 

£ = ye^aLe aR <f + h.C. 

=> e aR + 0° < — > e aL . 

(vi) The gauge interaction of the doublet Higgs fields : C = ^(5 M 0~)0°H / M h + h.c. 



(44) 



(45) 



+ 0° * — ► W~ 



(46) 



The following interactions (vii),(viii),(ix) do not exist in the standard model. 



(vii) The Yukawa interactions between the leptons and the triplet Higgs fields in Eq. (J3J) : 



1 J L v ' s ^2 



£° < — * v aL + upL, £ + < — ► v aL + ep L , £ ++ < — > e aL + e pL . 



+ h.c. 

(47) 



viii) The cubic coupling between the doublet Higgs and the triplet Higgs fields in Eq. (jSj) : 



Qcubic 



0°) 2 r - v^vr + (0 + ) 2 r- 
e — o" + 0°, e — 0° + + , e + — + + 



h.c. 



(48) 



(ix) The gauge interaction of the triplet Higgs fields : £ = g[(<9 M £ )£°— (<9 M £ )£ + }W++h.c. 

=> r + e — w~, c-^r + w-. (49) 

There also exist the sphaleron-anomaly processes through the non-perturbative effects 
in the standard model. 

(x) The sphaleron-anomaly processes 



Ui L u iL di L < — > e aL , u iL d iL d iL < — > u aL . 



(50) 
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We are ready to write down the Boltzmann equations in the SU(2)l triplet Higgs model. 
We first write down the Boltzmann equation for the time evolution of the lepton number 
L. The change of the lepton number L is caused by the lepton number violating processes, 
which are the Yukawa interactions between the leptons and the triplet Higgs fields (vii) and 
the sphaleron- anomaly processes (x). We write down the Boltzmann equation for the time 
evolution of the lepton number density n\ in Eq. (jA2|) under the assumption that the phase 
space density of a particle X is given by the Maxwell-Boltzmann distribution, 



In Eq. f)A2|) . Mf(f aa ) is the invariant matrix element for the processes, £° < — > v a i + v&l in 
(|47jl . The other invariant matrix elements, Mx are defined in the similar way. Although we 
can not exactly calculate the invariant matrix element for the sphaleron-anomaly processes 
M s , the obtained results in this section do not depend on the detail of the element M s , 
because we require that the sphaleron-anomaly processes are in equilibrium at large tem- 
perature region later. We can also write down the Boltzmann equation for the anti-lepton 
number density n\ in the similar way. The difference of the two equations leads to the the 
Boltzmann equation for the lepton number density n^[=ni — nj) 



f{E x ) = exp 



Ex — l^x " 
T . 



(51) 



n L + 3Hn L 



\M f {D\ 2 {-W - 2,0 + 2\M f (D\ 2 (-w ~ AS ~ As) 



+2|M / (/ er )| 2 (-/i ? o - ^ - /v) + 2\M f (D\ 2 (-W - AS - AS) 



+ 



\M f (f aa )\ 2 (^- ~ AS - ^l) + \M f {r)?W- - AS - AV) 



+ \M f (D\ 2 (^- ~ AS ~ AS) + \M f (D\\^- - ^ - /i T J 



+ \Mf(D\ 2 (^ 
+ \Mf(D\ 2 (^ 



AS - AS) + l^/(/ er )| 2 (^- - AS - AS) 



AS AS 



+ 



\M f (n)\ 2 (^- -2^ L ) + 2\M f {D\ 2 {^- -AS~AS) 



+2\M f {r)\ 2 W- -n eL -ti TL ) + 2\M f (D\ 2 (^- 
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du J A {PX - Pi - P2 - Ps) 



S 4 (p x ~Pl - pi) 

1 

2 



Y -ir\ M s\ 2 ( 2 ^u L +Hd L +V« L +Vu L + 2n dL + n Va ) , (52) 



where h L = dn L {t)/dt, f| = / rfn x rfn 1 rfn 2 (27r) 4 5 4 (px - Pi - p 2 ), fdTU = J fSpafe, and if 
is the Hubble parameter. Here we use the following relations, 

(53) 





■ E — fi X ~ 




" E + fix' 


E 


2/ijr 


exp 


I T . 


— exp 


I T . 




T 



/(l)/(2) - /(l)/(2) ~ e-f ■ ^i±M, (54) 

/(l)/(2)/(3) - /(l)/(2)/(3) ~ e-f ■ 2(/il+ ^ 2 + /i3) , (55) 

where the approximation fii/T <C 1 is used. 

We next write down the Boltzmann equation for the baryon number B. The change of the 

baryon number B is caused by the baryon number violating process, which is the sphaleron- 

anomaly processes (x) only in this model. We write down the Boltzmann equation for the 

baryon number density n& in Eq. (|A3|) . We can also write down the Boltzmann equation 

for the anti-baryon number density n^. The difference of the two equations leads to the 

Boltzmann equation for the baryon number density %(= - ng) 

f E 2 

h B + 3Hn B = / dnxdn 1 dU 2 dU 3 (27r) 4 5 4 (px - Pi - P2 - Pa) ■ e _7f • - 

Y -yl M s| 2 0"<*i. + 2 Vu L + ^d L + +Hu L +2fi dL ). (56) 

a=e,fi,r 

We write down the Boltzmann equation for the triplet Higgs number A. The change of 
A is caused by the interactions (vii) and (viii). We write down the Boltzmann equation for 
the triplet Higgs number density n$ in Eq. (|A4j) . We obtain the Boltzmann equation for 
"a(= n 5 - rig) 

h A + 3Hn A = p| -e-f • — 

[ Y ^|M / (r Q )| 2 (-2^ Q - / i C o) + |M / (r^)| 2 (-^ e -^-^o) 

+\M f (D\ 2 (-^ e - ii VT - Me o) + \M f {D?{-^ - /V - M?o) 



+ y \\MfU aa )\\-^ a -^L + ^) 



2 

a=e,fi,T 
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+ l\ Mf (D\ 2 (-»u e - + + \\M f {r)W-^ e - HrL + A*C") 

+\\M f {D\\-^ - HrL + + l -\M f {D\\-^ L - ^ + K-) 



+ Y \\M f U aa )\\-^ aL + ^..) + \M f {D\\-^L- n v 



+ \M f (D\ 2 (-VeL ~ HrL + ft-) + \M f (f 



P T \ I 2 / 



+ -\M A {A)\ 2 {2ii 4>0 - up) + -\M a (A)\ 2 (h4>° - H~ + ^-) 
+ l -\M A {A)\\-2 H -+^- 



(57) 



We similarly write down the Boltzmann equation for the doublet Higgs number <£>. The 
change of $ is caused by the interactions (iii), (iv), (v), and (viii). The charged current 
processes in companion with (iii), (iv), (v) also contribute to the changes of the number <3>. 
We write down the Boltzmann equation for the doublet Higgs number density in Eq. 
(IA5J) . The Boltzmann equation for n$(= — n$) is 

2 

n$ + 3Hn<f, = ("} -e~T ■ — 



Y \ M d{Vd)\ 2 {-^d lR + Hd lL - /ty>) 

i=d,s,b 

+ Y \ M u(y l u)\ 2 (vu m - Hu iL - /vO 

i=u,c,t 

+ Y \ M e(Ve)\ 2 {-V<xR + HaL- fijfl) 



i 

+ Y\ M u(Vu)\ 2 ^ - Vd iL + H<t>-) 
i 

+ y \M e (y:)\ 2 (-fi aR +^ a +^- 



a=e,fj,,r 



+\M A {A)\ 2 [{-^- + 2^-) + (-/x r - /i o + n<t>-) 

+ (/^o - 2//0o)] 



(58) 
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We finally write down the Boltzmann equation for the number of the charged gauge boson 
W. The change of W is caused by the interactions (i), (ii), (vi), and (ix). We write down 
the Boltzmann equation for the number density n w - in Eq. (jA6j) . The obtained Boltzmann 
equation for n^(= n w - — n w +) reads as 

E 2 

h w + 3Hn w = p| -e~T ■ - 



I Y \ M 9(A^ v km)\ 2 {-^u iL + Vd jL - Mu>-) 

i,3 



a>=e,fi.T 



+ \M g (—k^)\ 2 (^ <t> ~ 

+ |M 9 (5(A; M )| 2 ( / u 5 - + /i € o - 

+|M 9 (<?/g| 2 (/i e -- - // r - ^-)] . (59) 

Since the SU{2) L triplet Higgs model is within the framework of the SU{2) L ® U(l)y 
gauge theory, the total electromagnetic charge Q and the total third component of the 
SU(2)l isospin 1% are conserved. It is natural to assume that the total quantity of the 
conserved charge associated with the gauge symmetry vanish. Under this assumption Eqs. 
(1301) and (jH)) reduce to the following relations, 

Q = 0^ - 4/i € — - 2/i f - - A/i w - + 2 Y (VuiL + Vu iR ) 

i=u,c,t 

- Y (vd iL + vd iR ) - Y ^ aL + ^r) = ° ( 6 °) 



and 



+^( Y vu iL - Y vd iL ) + \ Y - = 0. (6i) 

i=u,c,t i=d,s,b a=e,fj,,T 

We can regard the interactions (i) ~ (vi) and (ix) as being in enough equilibrium. Since 
the creation and annihilation rates of a particle in in-equilibrium processes are equal, the 
in-equilibrium processes do not change the number of a particle. This observation yields the 
following relations about the chemical potentials from (|32*Jl . (jHH), (jSJ), (f15|). (jlSJ), and 



^u R = ^u L + /V 3 ) ^a R = ^u a + — ^<j)0, = [l m - — /i</,0, (62) 
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= -ftp, 




= -up, 





where we take the seven parameters, ((jl w -, \i Ul , \i Va , /x^o, \ip) as the independent parameters 
out of the nineteen ones. In Eq. (|62j). since the in-equilibrium processes (I42|) mix the up- 
type flavors of the left-handed quarks enough, the chemical potentials of them become all 
equal : fi UiL = fi UL (i = u,c,t). Furthermore, since the mixed flavors are transmitted to 
the up-type flavors of the right-handed quarks through the interactions (J44)) . the chemical 
potentials of them also become identical : \i UiR = \i UR (i = u, c, t). The similar processes lead 
to the similar relations for the chemical potentials of the down-type quarks as fid iL = ^d L 
and fid iR = /J-dR (i — d, s, b). Eq. (jpTj) restricted by Eqs. is transformed to the relation, 

Hu- = 0- (63) 
Inversely, Eqs. (J52~j) restricted by Eq. (jHSj) are written as, 



(64) 



where the six parameters fJL Ua , are taken to be the independent ones. Eq. (JfiUj) 

can be written by the six independent parameters as 

7/i0o + 3/ip + 6fi UL - 2fi u = 0, (65) 

where we define \i v as 

fi u = /W (66) 

Eq. (|63|) makes the number of the independent parameters reduce from six to five. The right 
hand sides of the Boltzmann equations (J52j) . (JoTlj) . (|57J). (|55J). and (JoTfy can be rewritten by 
the six parameters as in Eq. (|X7|I . (EU), iftgj) . (IMojl . and (JXTT}, where we use the relation, 

dL(t) dn L {t) 

s— = ^ r+ 3n L (t)H. (67) 
For simplicity, we further assume that the three lepton flavors, e, /x,r are identical, that is, 

j?ee f TT f^ T ~ f 

and fi Ve = fi„ = fi Ur =>- n v = 3(jL Ua . (68) 
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Since Eqs. (J68|) makes the number of the independent parameters reduce from five to three, 
we can take (L, B, A) to be the three independent parameters : 

1 





1 




1 


/V = 


k 


4 


Jfe 




1 






fl^O = 


k 








1 






= 










1 






= 




L 24 


+ 



24' 



i- £ - -La 

18 12 . 



The Boltzmann equations (jA7l 
parameters (L, B, A) as 



(69) 
(70) 
(71) 

Al . (72) 
18 12 J v ; 

9j), and (|M0|) can be also rewritten by the three 



7 



L 



1 



d 
dt 



B 



n 



2-e-f 
T-fc 



^ —7m/ — |m s 



4'' i / 12 

_ 17 

12 ^ s 



m., 



-3m/ + |m. 



B 



where 



m f = \M f (f)\ 2 , 
m A = \M A (A)\ 2 , 
in,. = I r/ll.i 



5\p X -Pi ~P2 ~p 3 ) 1, M ,2 



5 4 (px-Pi-p 2 ) 2 
Let us note that the determinant of the matrix in Eq. (|73j) . 

-7m/ - |m s |m/ - ^m s -3m/ + |m 



(73) 

(74) 
(75) 
(76) 



-lm a 
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51 



-m/ • m A • m s . 



(77) 



-|m/ + \m A |m/ - |m^ -|m/ - 
The in-equilibrium sphaleron processes (J5U|) at the temperature region between 100 GeV 
and 10 12 GeV lead to 



0. 



(78) 



which makes the number of the independent parameters reduce from three to two. We take 
(L, A) to be the two independent parameters and can reduce the Boltzmann equation (J73|) 
to 




n 



2-e-T 3 



-42m 



L 



T -k 17 y -21m f + \m A -f m f - f|m A J \ A 
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(79) 



We rewrite the Boltzmann equation (|79*jl into a simple form, 



s^-N = IM 2 N, (80) 
at 



where we define N, I and Mo as 



3 / -42m/ -fm f 
and M 2 = ^ J _ 2 f . (82) 
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-21m/ + fm A -f m/ - §m A 



Note that the determinant of the matrix, M 2 , is 



/ 3 \ 2 289 

det(M 2 ) = (-J ._. m/ -m A . (83) 

Here, the baryon number is also written by the linear combination of the two independent 
parameters, L and A as 

B = —-L + —A. (84) 
51 17 V ; 



B. Solution of the Boltzmann equation 

In this subsection, we solve the Boltzmann equation (|8()jl . To solve Eq. (|8())1. we need 
to know the JM 2 as a function of the time t or the temperature T. The factor J • nif is 
rewritten as 

where the decay rate r(£° — > P e P M ) is defined as 

r«» s ^ / gL^^^wv^ -» - P2 mu r ,86, 



and the triplet Higgs field £o has the energy E = ^/\p(\ 2 + M 2 . We define the thermal 
averaged decay rate of the decay process, £° — > v e v^ as 



r, s <r«« P . P ,)> = M^-^L! (87) 

J (27T)3 e T 
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and rewrite the factor I ■ rrif as 



I-m f = ^-nZ(T)-T f , (88) 



where 



d 3 p£ 
(2^ 



"SCO S / ^ e - T . (89) 



Similarly, we obtain the factor / • as 



/ " ™a = ^ • «m( t ) " Fa, (90) 
where is the thermal averaged decay rate : 

r A = 2(r(e° - 0V°)>- (9i) 

Using the thermal averaged decay rates and we rewrite Eq. f!80|) as 

-421/ -f^ \ / I 



d j 








eft 1 




1 T 3 





3"^ 4"/ 12" A 

To solve the Eq. (f9~2~j) we need to know the temperature dependence of the Tf and Ta- 
The decay rate r(£° — > v e Vp) is exactly calculated as 

r(£° ^ p p ) = ^! = £L (93) 

M 2-7T E E v ; 

We estimate in the two cases, (i)T > M and (ii)T < M. 

(i) T> M 

We do relativistic approximation, M ~ 0, which leads to 

J^r(e°-^ e )e-f _ aT2/27 r 2 



(94) 



= rfW) = T 3 /tt 2 = 2T' 

(ii) T < M 

We do non-relativistic approximation, ~ M + p 2 /2M and obtain 

/(0r(e°^^ e )e-f QW % /M o_ 
' „S(T) n% M' 1 J 
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The decay rate r(£° — > is exactly calculated as 



We similarly estimate in the two 

, /3/2T (T > M), 
T A ={ (97) 
/3/M (T < M). 

Since the temperature dependence of the right hand side of Eq. (j^j) is revealed, we 
rewrite the left hand side of the equation as a function of temperature T using the relation, 

s = "« T V (98) 



where we use the relation, 



Here the constant q is defined as 



H = qT 2 = i (99) 



q = lM^g:-^-, (100) 

Mpi 



with Mpi being the Plank mass and g* is the total degrees of freedom of effectively massless 
particles. For the following analysis, we introduce Kf and Ka as 



Kf =H 



and Ka = 1 



t=m 2qM 3 



(101) 



t=m 2qM 3 H 

Using the constants Kf and Ka-, we transform Eq. ()92|) into the Boltzmann equation which 
has the only one parameter T, 

-^N = f(T) ■ MN, (102) 
where the matrix M and the function f(T) are defined as 



2*/ T4 K f ^ 


l-l 


' a d \ 






{ cb) 



M " ' X A" "A"-- |S ' ' (1 ° 3) 
7 # (T > M) 



, M 3 

»»<i n-n { 1; don 

T 9/2 



^e-f (T < M) 
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The constants 7 and 5 in Eq. ()1U4|) are denned as 

378 1 , r 378 1 , , 

7 = — • — ~ 2.26 and 8= — - — = ~ 2.83. 105 

' 17 7T 2 17 7T 3 /2 .y/2 V ' 

We can solve Eq. ()102|) by diagonalizing the matrix M. The equation of the eigenvalue 
determines the two eigenvalues as 

\M-XI\ = Q x= a + b±D ^ (iog) 

where 



D = y/(a - b) 2 + Acd. (107) 

Using the transformation, 

/ X \ . ^ [2d -2d \ 

N = = V^N, with V = , (108) 

yY J yD + b-a D + a-b J 

Eq. ()102j) is transformed into the following one, 

±N' = f(T)-MN\ (109) 
where the matrix M is diagonalized as 

1 Ai , 

M = V~ l MV = . (110) 



A 



2 

Here the two eigenvalues Ai and A 2 are specified as 

= a + i + D = a+t-D 

2 2 v ; 

We can integrate Eq. (|109|) in the two cases, (i)T > M and (ii)T < M. 
(i) T > M 

The equation for the eigen-mode X in Eq. (J109)) is integrated as 

dX x M 3 . , 

if = ^7*-* (112) 

« X(T) = XtTJexpf-^M 3 ^-^)]. (113) 

It is natural for us to assume that M Tj, where Tj is the very high temperature 
in the very early universe, for example, Tj ~ 10 15 GeV. Then, the solution (|113|) is 
written as 

X(T — M) — X(Ti)e- l i Xl1 . (114) 
22 



(ii) T <M 

The equation for X is integrated as 

dX 



M 7 ' 2 

X(T f ~ 0) = X(Tx = M)e" 



J"9/2 



-3.19<5Ai 



(115) 

(116) 
(117) 



where we set Tj ~ 0, which is about the present temperature in the universe, and 
Ti = M. 

We combine the solution ()114|) with solution ()117|) at the temperature T = M and obtain 
the time evolution of X from Tj to Tf ~ as 



X(T) ~ 0) = X(Ti)e 



-rAi 



where 



1 

r = -7 + 3.195-9.77. 
3 

The equation for another eigen-mode Y in Eq. (jl()9)l is similarly integrated and Eq. 
is solved in a simple form, 

N' f = e-^Nl. 

When we go back to the original base N, the solution is given as 



(118) 
(119) 

jmnj) 

(120) 



N 



f 



-rAi 







L; 



v- 1 

e- rAa / \ A 



;i2i) 



where the evolution matrix S is concretely calculated to be 
S = Ve-^V' 1 

1 ( (D + a- b)e- rXl + (D - a + b)e~ rX2 2d(e~ rXl - e~ rX2 ) 

2D \ 2c(e~ rXl -e~ rX2 ) (D - a + b)e~ rXl + (D + a - b)e~ rX2 

Thus the final lepton and triplet Higgs numbers are obtained as 



;i22) 



J f 



A, 



2D 



2D 



[(D + a- 6) Li + 2dAi}e- rXl + \{D - a + 6)L< - 2rfAi]e" rA2 
\2cU + (D~a + 6)A 4 ]e- rAl + [-2cL; + (D + a - 6)A i ]e" rA2 



(123) 
(124) 
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Using Eq. 
Bf 



at the temperature Tf, we obtain the final baryon number, 



1 

2D 
+ 



28 4 
— (D + a-b) + — c)L,4 



51 17 v ; 



28 4 

(D -a + b) c)Li 

51 V ; 17 ' 



56 2 

— d+— (D + a-&))A< 



-rA2 



(125) 



It should be noted that the determinant and trace of the matrix M are both positive 

289 



det(M) = =^ • K f ■ K A = Ai • A 2 > 
67 5 

and Tr(M) = —K f + —K A = X, + A 2 > 0, 



(126) 
(127) 



28 J 36 

which means that the both of the two eigenvalues, Ai and A 2 , are positive. The positive 
eigenvalues tell us that the lepton number violating processes never generate the baryon 
number. Eq. ()126)) also tells us that if either of Kf or K A is zero, the determinant of M 
is vanishing, that is, A 2 =0. The vanishing of the second eigenvalue warrants the quantity, 
P = Y oc (bL — dA), which is an eigen-mode in Eq. (|108jh to be conserved and guarantees 
that the initial baryon number surely survives in proportion to the conserved U(l) charge 
P as 



or K A = det(M) = 



\ 1 = a + b, A 2 = 

B fm cx P = Y a (bL - dA). 



(128) 



In each of the two cases, (1)K A = or (2)Kf = 0, the conserved U(l) charge P concretely 
becomes the following, 



(1) K A = => P = L - 2A 

(2) K f = P = L 



(129) 



C. Condition to protect the primordial baryon number 



In this subsection, we analyze the solutions (|123|1 . p24)l . and p25jl of the Boltzmann 
equation and obtain the condition that the initial baryon number is not washed out. From 
the solutions ()123|) . ()124|) . and ()125j) we immediately find that the condition, Ai < 1 and 
A 2 < 1 approximately conserves the lepton, triplet Higgs, and baryon numbers as 



Ai < 1 and A 2 < 1 



J f 



ztLi, Af 



Ai, and Bf ~ B, L , 



(130) 
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where the two approximately conserved charges, P = L and A independently exist. However 
in order to protect the baryon number, two conserved charges are not needed. The only 
one conserved charge is enough to protect the baryon number. Hence the condition that the 
initial baryon number is not washed out is the following, 

Ai < 1 or A 2 < 1 A 2 < 1 <£> B f o:P=(bL-dA)cxY, (131) 

where we use the relation, Ai > A 2 in Eq. (jlll)) and the quantity P is approximately 
conserved. Using Kf and Ka, we rewrite the condition, A 2 < 1 as 

a + b — D 



A 2 < 1 <£> < 1 <£> Tr(M) - 2 < v/[Tr(M)] 2 -4det(M), (132) 

where det(M) and Tr(M) are defined in Eq. ()126|) and ()127|) . Inequality ()132j) is solved in 
the following two cases, (i)Tr(M) - 2 < or (ii)Tr(M) - 2 > 0. 

(i) Tr(M) - 2 < 

The condition, Tr(M) — 2 < is transformed to 

67 5 

28 K > + 36 /G " 2 < °- (133) 

(ii) Tr(M) - 2 > 

Inequality ()132|) is transformed to 

^ / + i A -2>0 and ( Kf -^).(K A -^-)<c u (134) 
28 i 36 V i 578/ V A 578 / ' y ' 

where the constant c\ is 

/882\ 2 / 335 289\ n n/lo . 1QC , 

ci = ~ 0.043. (135) 

1 V289/ V1008 882/ v ; 

Either ()133j) or p34j) is the solution of inequality ()132j) and is shown in Fig. ^ 
The condition A 2 < 1 is rewritten for the damping factor as 

e ~r\ 2 > e -9.77 _ 5 7 x 1Q -5 ( 13 g) 

If the initial baryon number is sufficiently large, the larger damping factor is allowed. For 
example, if we set the condition, A 2 < 2 to protect the baryon number, the damping factor 
is the following, 



e 



-rA 2 > e -2x9.77 _ 3 3 x 1Q -9 ^37) 
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FIG. 1: The region where the baryon number is not washed out is shown as a gray region. Here 
we require that the condition to protect the baryon number is A2 < 1. 

If the initial baron number is B>i ~ 1 as in the Affleck-Dine scenario, the hard damping ()137j) 
can realize that the present baryon number, B ~ 1CF 10 . The condition A2 < 2 is rewritten 

as 

Tr(M) - 4 < v/[Tr(M)] 2 -4det(M), (138) 

and the allowed region for Kf and Ka is shown in Fig. El It is appropriate that we set the 
condition to protect the baryon number to be 

K f < 1 or K A < 1, (139) 

for the various damping factor. When this condition is satisfied, the quantity P = bL — dA 
is approximately conserved. The relation between the condition ()139j) and the conserved 
number P = bL — dA is consistent with the relation (J128)) in the limit Kf — ► or K& — > 0. 
In the limit Kf —>■ or Ka — ► 0, the approximately conserved number P becomes the exact 
conserved number associated with the exact symmetry U(l)p. 

We next investigate how the final baryon number depends on Kf and Ka under the fixed 
initial conditions, Lj and Aj and confirm the validity of the condition ([139)1 . 
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FIG. 2: When we require that the condition to protect the baryon number is A2 < 2, the allowed 
region is shown as a gray region. 

We first observe the case that the initial condition is 

Li = 1 and A» = 0. (140) 

Under the initial condition (jl4U|) . the final baryon number ()125|) depends on Kf and Ka as 
in Fig. El It is found that the baryon number is not washed out in the region Kf < 1 or 
Ka < 1 in Fig. El We second observe the case that the initial condition is 

U = 1 and A, = 1. (141) 

Under the initial condition (j!41|) . the final baryon number ()125|) depends on Kf and Ka as 
in Fig. 01 In this case, we again find that the baryon number is not washed out in the region 
Kf < 1 or K A < 1 in Fig. H 

We finally study the case that the initial condition is 

Li = 1 and A { = -. (142) 

Under the initial condition (|142j) . the final baryon number (|125|) depends on Kf and Ka 
as in Fig. El In the figure, we find the anomalous behaver of the final baryon number. It 
is shown that the baryon number is washed out at the region Ka = 0, where the baryon 
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(Ki-Q.S4)(Ka- 14.6)= 0.17 



Region to protect 

the baryon number 




FIG. 3: The dependence of the final baryon number on Kf and Ka is shown when the initial 
condition Li = 1 and Aj = is selected. 

number should be protected. The reason is that the initial condition (J142)) corresponds to 
vanishing the quantity P = L — 2A which becomes a conserved charge at = and the 
final baryon number is proportional to P, that is vanishing. The similar behaver occurs 
when the initial condition, Lj = and Aj = 1 is selected. In this case, the baryon number 
is washed out at the region Kf = because the initial condition erases the quantity, P = L 
which becomes a conserved charge at Kf = 0. However, since the solutions ()123|) . ()124jh 
and (|125|) are not the exact ones of the Boltzmann equations ([52)1 . ()5b|). (|57j). (|58jl. and ()59p. 
the conserved charge P = bL — dA in (|128p is actually not the case. The true conserved 
charges are the following, 

K A = P = B - L + 2A, 

K f = P = B-L, (143) 

which are consistent with the conserved ones in the Boltzmann equation (|73*j) for the numbers, 
S, L, and A. The anomalous behaver found above actually should occur when a initial 
condition erases the true conserved charge P in ([143)1 . 
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FIG. 4: The dependence of the final baryon number on Kf and Ka is shown when the initial 
condition Lj = 1 and A* = 1 is selected. 

We here generalize the condition (j!39)l in order that we can apply the condition that the 
baryon number is not washed out to the various mass models of Majorana neutrino. The 
condition that the baryon number is not washed out in the mass models of Majorana neutrino 
under the in-equilibrium sphaleron processes is generally that if the effective lagrangian 
density is constructed only by the in-equilibrium interactions in the lagrangian density of a 
model, the effective lagrangian density has at least one conserved charge P associated with 
a group U(l)p symmetry which contains the lepton number. In other words, the condition 
is that when the out-of-equilibrium interactions are eliminated in the lagrangian density of 
a model, the lepton number is not explicitly broken in the truncated lagrangian density. 
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FIG. 5: The dependence of the final baryon number on Kf and Ka is shown when the initial 
condition Lj = 1 and Aj = 1/2 is selected. 

In the above discussions, the lepton flavor has been omitted because of Eq. ()68|) for 
simplicity. When we take the lepton flavors, e,/i, r into account, the condition that the 
baryon number is not washed out is extended. Since the sphaleron-anomaly processes exactly 
conserve the three independent charges, 

B B B , , 

y-^e, J~ L ^ aIld J~ L -> ( 144 ) 

the condition becomes that when the out-of-equilibrium interactions are eliminated in the 
lagrangian density of a model, at least one linear combination of the three charges (j 144)1 
is not explicitly broken in the truncated lagrangian density. Here the number of a lepton 
flavor, L a (a = e, //, r) is defined as 

iv a = 1 | 

L s s s 

Thus when we take the lepton flavors, e, /i, r into account and require the cases that the only 
one charge is approximately conserved, there are ten possible approximately conserved U(l) 

30 



(145) 



charges in the SU(2)l triplet Higgs model. We summarize all ten possible U(l) charges and 
the conditions to enhance the U(l) charges approximately in Table |l] In order to explain 
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TABLE I: In the SU(2)l triplet Higgs model, all ten possible U(l) charges and the conditions to 
enhance the approximate U(l) charges are shown. The condition K] Je < 1 means that when we 
write all the processes which violate the charge L e as Tl b , Tl £ < H\t=m- The other conditions 
have similar meanings. 

the condition K Le < 1 in Table HI we define T Le as the reaction rates of the processes 
where the charge L e is violated. For example, i\ e are the reaction rates of the processes 
as £° < — > u e L + VfiL and £ + < — > u e L + e + , which are all caused by the coupling constants 
f ea (a = e,fi,r). The condition K Le < 1 means that T Le < H\t=m- The other conditions 
have similar meanings. Here the charges L efl , L eT , and L^ T are defined as 

L e fi — L T — L e — L^, L eT = — L e — L T , and L^ T = L e — — L T . (146) 

The condition that the initial baryon number is not washed out is modified from the condition 
COHi into 

"at least one of the ten conditions in Table H] is satisfied in the model". (147) 
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IV. SU(2) L TRIPLET HIGGS MODEL COMPATIBLE WITH KAMLAND AND 
WMAP 

We obtain the conditions (|147|) that the primordial baryon number is not washed out 
in the SU(2)l triplet Higgs model. In this section, we investigate whether the constrained 
model by the conditions (|147|) simultaneously satisfies the recent results of the neutrino 
oscillation experiments and WMAP, and the constraints on the p parameter in the LEP. 
In subsection IIV Al we first summarize the present status of the neutrino oscillation experi- 
ments. In subsection IIVB| we examine the condition (0) Ka < 1 in ten conditions in Table 
01 In subsection IIV CI we thirdly examine the other nine conditions (1) ~ (9) in Table 
In subsection IIV D[ we finally summarize the results of the present section and show the 
allowed and forbidden regions of the parameters in the model. 



A. Present status of the results from the neutrino oscillation experiments 

We summarize the present status of the neutrino oscillation experiments to make clear 
what constraints the results of the experiments give to the neutrino mass matrix. The 
Super-Kamiokande Collaboration(S-K) shows that there exist a mass squared difference A a 
and a mixing angle 6 atm in order to explain the atmospheric neutrino oscillation 

1.6 x KT 3 eV 2 <A a < 4.0 x KT 3 eV 2 , 0.88 < sin 2 29 atm < 1.0 (90% C.L.), (148) 

with the best fit values A a = 2.5 x 10~ 3 eV 2 and sin 2 29 atm = 1.00. The global analysis of 
the first results of KamLAND combined with the existing data of solar neutrino experiments 
shows that there exist a mass squared difference A s and a mixing angle 8 Q to account for 



the solar neutrino oscillation 



s squared diffi 

DDDDQ 



5.1 x 10~ 5 eV 2 < A, < 9.7 x 10~ 5 eV 2 , 1.2 x 10~ 4 eV 2 < A s < 1.9 x 10" 4 eV 2 , 

0.29 < tan 2 6 Q < 0.86 (3a level), (149) 

with the best fit values, A s = 6.9 x 10~ 5 eV 2 and tan 2 9 Q = 0.46. The CHOOZ experiment 
has put the upper bounds on the mixing angle # 13 p, 

sin 2 2^i3 < 0.1. (150) 
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Using the Maki-Nakagawa-Sakata (MNS) matrix U, we can diagonalize a Majorana neu- 
trino mass matrix M. v to Ai u , 



( mi ^ 



My 



\ 



m 2 
m 3 



xFMuU, 



(151) 



^2 

where the weak eigenstates in are understood to be the partners of the mass eigenstates 
of the charged leptons and contain the neutrino mass eigenstates. Since we assume that 
is a real symmetric matrix, we can parameterize the MNS matrix U as a orthogonal 
matrix, 



U = 



'10 

C23 S23 
— S23 C23 
Rmi—QmiR 



\ 



( 



C13 Si 3 

1 

-«13 C13 



^ ( cn s 12 ^ 
-S12 C12 
1 



23j- f l-13l< 7 13y 



12; 



(152) 



where c,- 



cos 6ij and = sin^. Using these parameters, the two mass squared differ- 



ences and two mixing angles in Eq. (|148|) and (|149j) are written as 



2h 



A, 



rrir, 



m 



9, 



at i n 



9 



23- 



e Q = e 



12- 



(153) 



The experimental data uDDDO for the solar neutrinos also show that the sign of m 2 , — m\ 
is positive. For brevity, we assume the mixing angle #13 = 0, consistent with ()15U|) . and fix 
the two mass squared differences and 9 a tm at the best fit values. Under these assumptions, 
we can write the matrix U using the only one relatively poorly known parameter 9 & = 9, 



( 



U 



cos 9 



sin 9 







\ 



-7I sin ^ 75 cos ^ 71 
V 75 sin ^ "7I cos ^ 73 / 



with 0.29 < tan 2 9 < 0.86 (3a level). (154) 



Then, the neutrino mass matrix M. v can be written as 

^ c 2 mi + s 2 m 2 — ^jsc(mi — m 2 ) -^sc(m\ — m 2 ) 

■M u — — -i=sc(mi — m 2 ) \(s 2 mi + c 2 m 2 + m 3 ) — |(s 2 mi + c 2 m 2 — m 3 ) 
, -^sc(m! - m 2 ) -|(s 2 mi + c 2 m 2 - m 3 ) |(s 2 mi + c 2 m 2 + m 3 ) , 



(155) 
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where c = cos# and s = sin#. Since we have fixed the values of the two mass squared 
differences, there remains only one independent parameter among three mass eigenvalues, 
m 1} m 2 ,m 3 . 

While the neutrino oscillation experiments do not determine the magnitudes themselves of 
the neutrino masses as in Eq. (|153|) . the investigation into the cosmic microwave background 
radiation (CMBR) byWMAP recently put the following upper bound on the sum of the 



neutrino masses 



10 



11 



3, 0,0, 



Q u h 2 = 54^T7 < 0.0076 <£> V \mA < 0.71 eV. 
93.5 eV ^-^ 



(156) 



In the following analyses, we require that the constraints (j!55J) and ()156p are satisfied in the 
model. 



B. First condition Ka < 1 

In this subsection, we examine whether the condition (0) Ka < 1 of ten conditions is 
compatible with the results (|155J) of the neutrino oscillation experiments, the results ()156|) 
of WMAP, and the constraints (J29j) on the p parameter. The condition (0) Ka < 1 in Table 
Uis rewritten from Eq. (jlOlj) as 

\A\<12x (^) 5 eV. (157) 
The constraints (|29j) on the p parameter is written as 

\ M < 0.03 x (158) 

V 



We compare the right hand side of ()157|) with one of ()158j) as 

( M \5 , r %V2M 2 
,GeV, 



12 x — - eV < 0.03 x 
VGeV/ 



<(=> M > 1.2 eV, (159) 

which shows that ()157j) is more severe condition than ()158|) because of M ^> 1.2 eV. We 
can easily show that the condition ()157|) is compatible with the results ()155j) of the neutrino 
oscillation experiments. When we fix the mass M at 100 GeV, for example, the condition 
(I157|) becomes 

A < 12 KeV. (160) 
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When we further assume that the mass spectrum is normal type and mi = 0,rri2 > 0, and 
m3 > and fix the solar mixing angle, 9 in (j!49|) at the best fit value, the mass matrix (jl55j) 
becomes 



/ 



Mi 



A 



( 2.6 x KT 3 2.7 x 10~ 3 -2.7 x KT 3 ^ 
2.7 x KT 3 2.5 x 10~ 2 2.5 x KT 2 
-2.7 x 10~ 3 2.5 x 10~ 2 2.5 x KT 2 



eV. 



(161) 



When we fix A at 1 eV, which is consistent with the condition ()160|) . f af3 is calculated from 
the equation, m a/3 ~ 1.5 eV x/ a/3 in (j2T?j) and ()161f) as 



./ 



a/3 



1.5 eV 



V 



1.7 x 1(T 3 1.8 x 1(T 3 -1.8 x 10 

1.8 x 1(T 3 1.7xl0~ 2 1.7 x 1(T 2 
-1.8 x 10~ 3 1.7 x 1(T 2 1.7 x 10~ 2 



-3\ 



(162) 



Inversely, when we take f af3 as (|162j) . both of the mass matrix (J161)) and the condition 
Ka < 1 are simultaneously satisfied. It is also easily found that the mass matrix (J161|) 
satisfies the result (|156|) of the WMAP because of the equations, 



\mi\ = = TrM v ~ 5.3 x 1(T 2 eV < 0.71 eV. 

i i 

In the case < 1, the approximately conserved charge is from (J143j) . 



(163) 



P = B - L + 2A 



(164) 



and the final baryon and lepton numbers are calculated from a equation, 

/i 5 o + 2fi„ a = 0, 

(J7UJ), (JI2J), and (jHU) and are protected in proportion to the initial value Pi n i as 



(165) 



B 



fin 



--L 
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We next investigate the region of the coupling constants f a/3 where the results of the 
neutrino oscillation experiments are not satisfied under the condition < 1. We use the 



p. 



(166) 
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boundary values of the other nine conditions (1) ~ (9) in Table U] as the standards which 
classify the regions of the coupling constants f a ". We show below that the constrained 
neutrino mass matrix (|155|) can't satisfy any one of the nine conditions in addition to the 
condition Ka < 1. The nine conditions (1) ~ (9) are written as a form, 

M y 
.GeVJ ' 



If 



a/3 1 



< 4.3 x 10" 



x 



(167) 



The condition ()157|) and ()167j) derive the constraint on the magnitudes of the elements of 
the neutrino mass matrix in Eq. (|26J) as 



\m 



a/3 1 



I./" 



" 3| 4S<7.9x 10- 4 eV (=S). 
AM 2 v ' 



(168) 



If the mass matrixes are restricted by the inequality ()168|) for all a = e, fi, r, the matrixes 
are not able to induce the atmospheric mass squared difference, \fK^. = 5.0 x 10~ 2 eV in 
(I148j) . We next impose the condition ()157j) and the one of the left nine conditions in Table 
U] We enumerate the type of the mass matrixes M u which are associated with the exact nine 
symmetries which the approximate nine symmetries go to in the limit, — > as 



P? = L, 


x x 
x x 



Pi = h 
x x 

x x 



P{ = L T 
^ x x 0^ 
x x 




(169) 



pa 
r 2 



0x0 
x 
x 



P 



pc 



x 
0x0 
x 



x 
x 
0x0 



170) 



pa 
r 3 



M„ 



( 



\ 



L e + — L 7 



x 
x 

x x 



P b 



L e ~ + L T P3 — — L e + + L 7 



(171) 



V 



0x0 
x x 
0x0 



/ 



V 



x x 
x 
x 



\ 



/ 



Thus we examine whether the mass matrixes associated with the nine approximate symme- 
tries can satisfy the results of the neutrino oscillation experiments or not. We first examine 
the case that the mass matrixes have the normal mass spectra. 
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(i) Normal mass spectrum : \mi\ < \m%\ < \m%\ 

Choosing mi as a parameter, we can express m 2 and m% as a function of mi, 



m 3 = ±,/A a + A s + m 2 , (172) 



m 2 = ±^/A s + m 2 , (173) 
mi > 0, (174) 

where mi can be taken to be non-negative with no loss of generality thanks to the 
freedom of the re-phasing, z7£V = ±zz7£\ Every element of the neutrino mass matrix 
can be written as 



m ee = c 2 m 1 ±s 2 \J A s + m\, (175) 

-^=sc(mi =f yjA s + mfj , (176) 

- l - [s 2 m x ±c 2 sjA s + m\T\J&a + A s + m 2 ^ (177) 

- (s 2 mi ± c 2 ^A s + m\ ± ^ A a + A s + m 2 ) . (178) 



m eM = -m eT 



m" T 



2 

For the approximate symmetries, (P", P 2 a , P 2 b , P3 , P 3 b , P|), \m ee \ in Eq. ()175|) is re- 
stricted as 



|m ee | 



< 5 = 7.9 x 10~ 4 eV. (179) 



c mi ± s y A s + mf 

If m 2 > 0, Eq. (j!75J) induces the inequality, 

m ee > s 2 ^/A~ s > s 2 min ^A~ s = 1.8 x 10~ 3 eV > S. (180) 

Inequalities ()179j) and ()180)) need the inequality m 2 < 0. Here we shorten sin# m j n as 
s min and 6 min is defined as tan 2 6 min = 0.29 in Eq. ffHHj). For (Pf, P 2 a , P 2 6 , P 3 a , P|), the 
element |m eAt | or |m eT | has to be restricted as 

|m e ^| < S or |m er | < 5. (181) 

On the other hand, Eq. ()176j) with m 2 < derives the inequality, 

|m eM | = |m eT | = -^=sc^mi + A s + mfj 

> -=s min c mi „^A; = 2.4 x 10~ 3 eV > S, (182) 
V2 
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where c m i n = cos Q m %n- Since f!182|) is inconsistent with ()181|) . the approximate symme- 
tries, (Pi, P£, P|, P£, P|) are not allowed. For P 3 C , the inequality |m MM | < 5 is needed. 
On the other hand, Eq. (|178|) derives the inequalities, 

m 3 > «=> m w > 2.2 x 10" 2 eV > S, (183) 
m 3 < <=> m w < -2.8 x 10" 2 eV < -5. (184) 

which don't allow the inequality \m m \ < S, that is, the approximate symmetry Pg. 
The left approximate symmetries, (Pf, Pf, P 2 C ) need the inequality, 

|m e ^| < 5 or |m er | < 5, (185) 

which require rri2 > because of ()182|) . The symmetries further need the inequality, 

\m^\ < S or |m rr | < S. (186) 

On the other hand, Eq. ()178|) with m 2 > derives the inequality, 

m 3 > => |m wt | = |m rT | > > S. (187) 

Since the inequalities ()186|) and ()187j) are inconsistent, the inequality m 3 < is needed. 
For Pf and Pf, the inequality, \m^ T \ < S is needed. On the other hand, Eq. ()177|) 
under the inequalities m 2 > and m 3 < induce the inequality, 

\m^\>^^/A~ a > S, (188) 

which don't allow \m^ T \ < S, that is, the approximate symmetries Pf and Pf. The left 
one approximate symmetry P| needs the inequality, |m w | < 5 which is solved under 
m 2 > and m 3 < for \mi\ in Eq. f|178|) as 

\ m w\ < 5 =>• |mi| > 0.80 eV, (189) 

which is inconsistent with the results (I15fi|) of WMAP and don't allow the approximate 
symmetry P|. From the inconsistency obtained above, we find that none of the nine 
approximate symmetries can be compatible with the neutrino mass matrixes with 
normal mass spectra under the condition Ka < 1. 
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(ii) Inverted mass spectrum : |m 3 | < |mi| < \rri2\ 

Choosing the smallest mass m 3 as a parameter, we can express mi and m 2 as a function 
of m 3 , 



m 2 = ±^A a + A s + m 3 , (190) 

mi = ±y/A a + m 3 , (191) 

m 3 > 0, (192) 

where m 3 can be taken to be non-negative without loss of generality. Every element 
of the neutrino mass matrix can be written as 



m ee = ±c z s JA a + m| ± sy A a + A s + m|, (193) 
= -mT = __L sc (±^ Aa + m 2 T ^A a + A s +mt), (194) 



-- (±s 2 ^A a + ml±c 2 ^A a + A s + ml- m 3 ) , (195) 
l - (±s 2 sjA a + m\ ± <?yj A a + A s + ml + m 3 ) . (196) 



For the approximate symmetries (P?, P£, P%, P 3 , P 3 , P 3 ), the inequality |m ee | < S is 
needed. On the other hand, |m ee | under the inequality, mi ■ m 2 > in Eq. ()193j) 
derives the inequality, 

|m ee | > v^> S, (197) 

which is inconsistent with |m ee | < 5 and require the inequality, mi • m 2 < 0. In the 
approximate symmetries (P", P 2 a , -^t)' the element |m eM | or |m er | has to satisfy 

the inequality, 

|m^| < 5 or |m er | < S. (198) 

On the other hand, Eq. ()194|) with the inequality, mi ■ m 2 < derive the inequality, 

\ m ^\ = | m -| > V2s min c min ^fA~ a = 2.9 x 1(T 2 eV > S, (199) 

which is inconsistent with the inequality ()198j) and doesn't allow the approximate 
symmetries (P? , P£ , P% , P^ , P^) . For P|, the inequalities, |m w | < S and (m^l < S 
are needed. On the other hand, Eqs. ()195|) and ()196|) derive the inequalities, 

mi > and m 2 < ^ |m Mr | > 2.0 x 1(T 2 eV > S, (200) 
m 1 < and m 2 > \m^\ > 2.0 x 10~ 2 eV > S, (201) 
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which are inconsistent with |m w | < S and \m tiT \ < S and don't allow the approximate 
symmetries P| . The left symmetries (Pf, Pf , P 2 C ) need the inequalities, 

|m e ^| < 5 or \m eT \ < S, (202) 

which require the inequality, mi ■ m 2 > from the inequality (|199|) . The symmetries 
Pf and Pf also need the inequalities, 

"|mH < 5 and \m^\ < S" 
or 

"|m Mr |<5 and |m Tr | < 5". (203) 

On the other hand, Eq. ()196|) and ()195|) derive the inequalities, 

1 , — 

mi > and m 2 > ^ \m^\ = \m TT \ > -VA a > 5, (204) 

mi < and m 2 < |m^ r | > > S, (205) 

which are inconsistent with the inequalities ()203|) and don't allow the approximate 
symmetries Pf and Pf. The left one approximate symmetry P| needs the inequalities, 
\m^\ < S and \m TT \ < S, which derives the inequality, vii\ < and m 2 < from 
(|2()4jl . The inequality, |m MM | < S with mi < and m 2 < is solved for from (J196j) 

as 

| TO w| < 5 ^ | m3 | > 0.81 eV, (206) 

which is inconsistent with the results f!156|) and don't allow the approximate symmetry 
P^. From the inconsistency above, we find that none of the nine approximate symme- 
tries can be compatible with the mass matrixes with inverted mass spectra under the 
condition Ka < 1. 

The above analyses show that both of the condition < 1 and any one of the nine con- 
ditions (1) ~ (9) in TableHJare not simultaneously satisfied in the mass matrixes compatible 
with the results of the neutrino oscillation experiments. Such results favor the relation that 
the coupling constant A is much smaller than the electroweak scale but the coupling con- 
stants f a/3 are the same order as the Yukawa couplings of the charged leptons in (J2*7|) . Such 
allowed region is expected to be natural in the SU(2)l triplet Higgs model in the framework 
of the large extra dimension. 
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C. Other nine conditions 



In this subsection, we impose the other nine conditions (1) ~ (9) in Table U] as 



(207) 



Here, what a and j3 satisfy the inequality ()207|) depends on the condition which we select 
in the nine conditions. Since in the previous subsection we find that both of the condition 
Ka < 1 and any one of the nine conditions are not simultaneously satisfied in the mass 
matrixes compatible with the neutrino oscillation experiments, under the condition ()207|) 
we need to search for the parameter region Ka > 1, that is, 



w > 12 * (^) eV - 



(208) 



The constraints (j!58|) on the p parameter and the region 
12x 



are written as 
eV. 



{ JLy eV<lA]< , 5xwH JLf 

The inequalities ()207|) and (J158J) give the upper limit on the neutrino masses as 



(209) 



\m 



a/3 1 



1/ 



a/3 



\A\v 2 
AM 2 



< 26x 



M \\ 



GeV 



) 



eV, 



(210) 



which seems to be compatible with the neutrino oscillation experiments. We actually show 
that both of the conditions, (J207j) and (|209|) are compatible with the results of the neutrino 
oscillation experiments in a typical example. When we assume that the neutrino mass matrix 
is (jl61J) and fix the mass M at 100 GeV again, the conditions ()207|) and ()209|) are written 



as 



\f af3 \ < 4.3 x 10~ 8 , 

8.6 KeV < \A\ < 3.5 GeV, 



(211) 
(212) 



When we further fix the coupling constant A at 100 MeV, the coupling constants f a/3 are 
obtained from the equation, m al3 ~ 151 MeV x/ a/3 and (J161)) as 



./ 



a/3 



151 MeV 



1.7 x 10- 11 1.8 x 10- 11 -1.8 x 10 

1.8 x 10~ n 1.7 x 10~ 10 1.7 x 10~ 10 
-1.8 x 10~ n 1.7 x 10~ 10 1.7 x 10~ 10 



-ii \ 



(213) 
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which are consistent with ()211|) for all a and (3. Hence all the nine conditions (1) ~ (9) 
are simultaneously compatible with the results of the neutrino oscillation experiments. In 
the case that all the nine conditions are simultaneously satisfied, the allowed approximately 
conserved charge is from (|143|) . 

P = B-L (214) 
and the final baryon and lepton numbers are calculated from a equation, 

2/i^o - 2/i € o = 0, (215) 
(jSHJ), (fTUj) . and (|S4^) and are protected in proportion to the initial value Pi n i as 

Bfin = ~J^Lfin = J^Pini- (216) 

Since we expect that the small coupling A and the Yukawa couplings f af3 , being same 
order as those of the charged leptons, to be natural in this model, the conditions (J211j) and 
(I212|) are necessarily regarded as the unnatural parameter region in the model. 



D. Allowed region compatible with all requirements 

In Subsection IIVBI and IIV CI we examine whether each of ten conditions in Table H] is 
compatible with the results of the neutrino oscillation experiments and WMAP, the con- 
straints on the p parameter and we obtain the allowed region of the parameter. Since the 
allowed region complicatedly depends on the six coupling constants f a/3 , the allowed region 
are not completely obtained. Here, for simplicity, we consider the only two cases in which 
all the coupling constants f a @ satisfy the condition ()2()7)1 and in which none of f af3 satisfy 
(120 7|) . We further fix the masses of the triplet Higgs fields at M =100GeV. Under the as- 
sumptions we show a typical allowed region in Fig. |U1 The three forbidden regions where the 
requirements are not satisfied are shown in the gray region in Fig. |BJ The conditions that 
the primordial baryon number is not washed out are not satisfied in the forbidden region (i) 
where both of the coupling constants A and f al3 are large. The results of the neutrino oscil- 
lation experiments are not satisfied in the forbidden region (ii) where both of the coupling 
constants A and f af3 are small. The constraint on the p parameter is not satisfied in the 
forbidden region (iii) where A is large. Although we exactly obtain the forbidden region, 
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FIG. 6: Three forbidden regions (i), (ii), and (iii) are shown as gray regions. Two allowed regions 
(1) and (2) which we confirm are shown as white regions. The masses of the triplet Higgs fields 
are fixed at M =100GeV. 

we do not exactly obtain the allowed region. What we show is that there exist the allowed 
region in the two white regions in Fig. The allowed region (1) in the white region where 
the coupling constant A is the very small and the Yukawa couplings / Q/3 are as large as those 
of the charged leptons in (|27j) is natural parameter region of the SU (2) l triplet Higgs model 
in the large extra dimension. Inversely the allowed region (2) in the white region where A 
is the large and f a " are very small is expected to be unnatural. 
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V. SUMMARY 



A. Conclusions 

We introduce the SU(2)l triplet Higgs model where the smallness of the neutrino masses 
is naturally explained by the shining mechanism in the large extra dimension. In order to 
violate the lepton number explicitly, which is the source of the Majorana mass terms, we need 
to introduce the two kind of the interactions in addition to the ones in the standard model. 
It is a first fruitful consequence that we can obtain the condition (|147|) that the primordial 
baryon number is not washed out by solving the Boltzmann equations in the model. We 
find that in order to protect the primordial baryon number, either of the two kind of the 
interactions, both of which are needed to violate the lepton number explicitly, must be out of 
equilibrium under the in-equilibrium sphaleron- anomaly processes. The condition to protect 
the primordial baryon number is easily extended into the case that the differences of the 
three lepton flavors are taken into account. It is a second fruitful consequence that we can 
obtain the allowed region where the condition to protect the primordial baryon number, the 
results of the neutrino oscillation experiments and WMAP, and the constraints on the p 
parameter are all satisfied in the model. A typical allowed region is shown in Fig. |U1 We 
find that if both of the coupling constant A and f a/3 are large, the existing baryon number is 
washed out and when both of A and f af3 are small, the mass differences which are induced 
from the neutrino mass matrix are too small to explain the observed neutrino oscillations. 
In this model, the constraints on the p parameter give the upper bounds on the coupling 
constant A for the fixed triplet Higgs masses. We find the two allowed regions where all 
the requirements are satisfied. One region where A is very small and f af3 are as large as 
the Yukawa coupling constants of the charged leptons is expected to be a favorable region 
where the smallness of the neutrino masses is naturally explained in the SU {2) L triplet Higgs 
model in the large extra dimension. 

B. Discussions 

We first make comment on the pre-existing approximately conserved charge. We actually 
obtain the approximately conserved charges, P = B — L + 2A or P = B — L. The primordial 
baryon number is protected in proportion to the initial value of the conserved charge. Can 
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we naturally expect that the initial amount of the conserved charge exists in the early 
universe ? We recall that in the SU(5) GUT, the initial value of the conserved charge 
(B — L) is naturally assumed to be zero, which regrettably leads to the washing out of the 
baryon number by the in-equilibrium sphaleron processes. The difference between the case 
we examine and the case of SU(5) GUT is that in the former case the conserved charge is 
approximate one, while in the latter case the conserved charge is exact one. So there is no 
reason why an approximately conserved charge is not generated by the unknown mechanism 
at the very early universe. 

We finally enumerate the remaining problems and the analyses which have not been 
carried out yet in this paper. The first problem is that we have omitted all four body 
processes such as 11 <-> 11 or 11 <-> 00, where I is the lepton and is the Higgs field. Although 
we expect that the rates of the four body processes are much smaller than those of the three 
body processes, the four body processes should be taken into account for the exactness. The 
second problem is that the analyses containing CP violation are not carried out. We have 
dealt in the restricted parameter region where there are no CP violation in the lepton sector 
of the models in order to focus on the feature that the baryon number is washed out in 
the Majorana neutrino mass models. If we try to investigate the feature more generally and 
more exactly, the analyses containing the CP violation should be carried out. However, even 
if the SU (2) L triplet Higgs model which has the only one triplet Higgs field contains CP 
violation in the lepton sector, the baryon number can not be generated. Hence the condition 
to protect the primordial baryon number which we obtain in this paper is valid in the case 
that the model contains the CP violation. On the other hand, since Majorana neutrino mass 
matrix generally has one Kobayashi-Maskawa phase and two Majorana phases, the analyses 
of the results of the neutrino oscillation experiments might be slightly changed. 
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APPENDIX A: BOLTZMANN EQUATIONS IN SEC. UTT] 

The difference between the number density of the matter and one of the anti-matter is 
calculated as 



n + — ?2_ = g 



d 3 p 



1 1 



e t ±\ e t ± l 



1 ,1 (fermion) 
= -gfiT x I (Al) 
6 ^2 (boson). 

We write down the Boltzmann equation for the time evolution of the lepton number 

density rii as 

m + 3Hm = J dn x dn 1 du 2 (2n) 4 5\ Px - Pl - P2 ) 

[ E l M /(/^| 2 [/(^)-/K)/K)]+2|MKni1/(f)-/(^e)/(^)] 
+2|M / (/-)l W) - / W/K)] + 2|M / (/^)| 2 [/(e°) - 



+ E l M /(/ aQ )| 2 [/(D-/K)/(e Q L)] 

+|M/(m 2 [/(r) - + iM/cm^cn - /(^/fo)] 

+ |M / (r T )| 2 [/(D - /(^/(Ti)] + iM^/^H/r) - /(e L )/K)] 
+ |M / (/ er )| 2 [/(D - f{e L )fW] + \M f (D\V(0 - /(Afc)/K)] 

+ E i^/(r a )H/(r-)-/(e Q L)/(e ai )] 

+2|M / (m 2 [/(r-) - /(ei)/(A*L)] + 2|M / (m 2 [/(r-) - /(e*)/fo)] 

+2|M / (/m 2 [/(r-) - /w/K)] 

+ / dll /fr*-***-^ (A2) 

E ^|M s | 2 [/(^)/(^)/(4) - /(e Qi ) + f(u L )f(d L )f(d L ) - f(v a )} " . 

We write down the Boltzmann equation for the time evolution of the baryon number 
density n b as 

h b + 3Hn b = J rfn x dn 1 dn 2 dn 3 (27r) 4 5 4 (p x - Pl - P2 - p 3 ) (A3) 
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E \\M s \ 2 U\e aL ) - f(u L )f(u L )f(d L ) + /(P Q ) - f(u L )f(d L )f(d L )]. 



2 

a=e,^t,T 



We write down the Boltzmann equation for the time evolution of the triplet Higgs number 
density ns as 

h s + 3Hn s = P| 

a=e,fi,r 

+\Mf(n\ 2 im)f(vr) - + w^mmm - f{e)\ 



+ E \\Mf(f aa )\ 2 lf(v a )f(e«L)-m + )} 

a=e,/i,T 

+\\M f {r)\ 2 \f{ve)f{H) - m + )\ + l\M f (n\ 2 im)f(f L ) - /(ni 
+l\M f (r)\ 2 [f(^)f(r L ) - + \\M f {D\ 2 \f^e)f{H) - f{e)\ 
+\\M f (n\ 2 [f^e)f(f L ) - /(ni + l\M f (D\ 2 m*MTL) - m + )\ 



+ E l\ M Af aa )\ 2 [f(e aL )f(e aL )-f(e + )] 

+ \M f (D\ 2 [f(e L )fm L ) - + \M f (r)\ 2 [f(e L )f(f L ) - /(£++)] 

+\M f {r)\ 2 [f{H)f{f L )-f{e + )] 



+ ^im a (a)| 2 [/(0°)/(0°) - f(e)} + i|M A (A)H/(0°)/(0 + ) - f(t)] 

+\\M A {A)\ 2 [f{ ( t> + )f{ ( t> + )-f{t + )]\ (A4) 

We write down the Boltzmann equation for the time evolution of the doublet Higgs 
number density as 

[ E \MM)\ 2 [f(d iR )f(d iL )-f(<P )] 

i=d,s,b 

+ E |m u (^)| 2 [/(M/(^l)-/(0 )] 

i=u,c,t 

+ E \Me(y a e )\ 2 [f\e aR )f(e aL )-f^)} 
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+ J2\ M M)W(d iR )f(u iL ) - f(<P + )] 

i 

+ Y^\ M M)\ 2 [f(u iR )f(d iL )-f(cj> + )} 

i 

+ E l M e«)| 2 [/(e^)/K L )-/(0 + )] 



+\M A {A)\ 2 \f{e)-f{<f)f{<f)} 
+\M A (A)ma-f(<f>°)f(<P + )] 

+ \M A (A)\ 2 {f(t + )-f(<P + )f(<P + )] ]■ (A5) 

We write down the Boltzmann equation for the time evolution of the charged gauge boson 
number density 

n w - + 3Hn w - = Q 

E \ A U^M)\V(^L)f(d jL ) - /(«,")] 

+ E \M g (j=)\ 2 [f(u a )f(e aL ) -f(w-)} 



+ \M 9 (gkXif(Om°)-f(w-)] 

+ \M 9 (gk ll )\ 2 lf(C~)f(t)-f(w-)}]. (A6) 

The existences of the in-equilibrium processes make the Boltzmann equations, (|52jl . 1)56)). 

(I57jl . (jSHI ) and (fo^j) the simplified ones as 

c?L __e 2 

s-r- = -e t ■ — 
dt 1 1 T 

[ E |M / (/n| 2 (-3)x( / i ? o + 2^J-6|M / (r^)| 2 (^o+ / i, e +^) 



5 4 (px - pi -Pz) 

E ^|M s | 2 (-2)x(3 /UuL+/ i QL )j, (A7) 



2 

a=e,ii,T 

dB f 2 

s— = / rfn x rfn 1 dn 2 rfn 3 (27r) 4 5 4 (p x - Pl - P2 - P3 ) • • - 
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a=e,/i,r 



(A8) 



_B 2 

■e T ■ — 
T 



[ E iM/(r)i 2 (y)(^+K)-3|M J (nf(^ +/ i,+^) 



mADW^e + + /O - 3|M / (/'* r )| 2 (^ C o + ^ + /i, T ) 



+ -|M A (A)| 2 (2^o -/x s0 ; 



(A9) 



_E 2 

e t . _ 
T 



-3|M A (A)| 2 (2^o-/i 5 o), 



(A10) 



W = 0. 



(All) 
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